UNIVERSALITY OF CORRELATIONS FOR RANDOM 
ANALYTIC FUNCTIONS 

SHANNON STARR 

^— ^ Abstract. We review a result obtained with Andrew Ledoan and Marco 

f"^ Merkli. Consider a random analytic function f{z) = y^^n anX„z'^, where 

^vj the X„'s are i.i.d., complex valued random variables with mean zero and unit 

^_^ variance, and the coefficients a„ are non-random and chosen so that the vari- 

r^ ance transforms covariantly under conformal transformations of the domain. 

h— -5 If the X„'s are Gaussian, this is called a Gaussian analytic function (GAF). 

We prove that, even if the coefficients are not Gaussian, the zero set converges 
in distribution to that of a GAF near the boundary of the domain. 



o 



Ph 1. Introduction 

''pi Random polynomials and random analytic functions are a topic of classical in- 

Cu terest. They have enjoyed renewed interest in the past decade or so. One can 

G imagine various motivations. The characteristic polynomial of a random matrix is 

'~^ a random polynomial. Similarly, the partition function for a Bernoulli spin glass 

, may be viewed as a random analytic function in the variable z = e^ . This is related 

►;>. to the celebrated Lee- Yang theorem. 

l/~j However, in mathematical study, it is easier to consider a different type of random 

C^ analytic functions, where the coefficients are independent. When the coefficients 

' ■ are i.i.d., the expected number of zeros has been well-studied. Consider for example 

. the random analytic function 
V 00 

ILJ where the coefficients are i.i.d. Then under general conditions, the zero set accu- 

mulates at the unit circle. A recent result |4j has found the sharp condition for 
the zero set to be asymptotically uniformly distributed on the circle. This type of 
result is akin to a law of large number: for example Ibragimov and Zaporozhets 
prove that the empirical distribution of the zeros converges, almost surely, to the 
uniform distribution on the circle. 

Various researchers have also considered the correlations between the zeros, un- 
der the assumption that the coefficients are all Gaussian. For instance, a recent 
monograph of Hough, Krishnapur, Peres and Virag reviews this topic [3]. A natu- 
ral question arises: if the zeros of a random analytic function accumulate near the 
boundary, even under the assumption that the coefficients are not Gaussian, what 
can we say about the correlations of the zeros near the boundary? For general rea- 
sons, one expects the distribution of the zeros to converge back to the distribution 
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of the zeros of the Gaussian analytic function, when properly rescaled. That is 
what we prove. 

One would expect to have various applications of such a result. However, some 
desired applications are presenly out of reach. The polynomials corresponding to 
random matrices and spin glasses generally do not have i.i.d. coefficients, except 
for companion matrices. There has been some recent interest in polynomials with 
discrete random coefficients. We will review this topic in the last section. 

Acknowledgments 
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2. Set-up 

We consider a special ensemble of random analytic functions, inspired by the 
monograph [3]. Given a parameter k < 0, and a sequence of coefficients x = 
(xo, xi, a;2, . . . ), one may define the power series 



where 



We consider random analytic functions (RAF's) defined by choosing a coefficient 
sequence X ~ {X^^Xi, . . .) where ATq, ATi, . . . are i.i.d., complex valued random 
variables, with mean zero and unit variance, such that 

(2.f) E[(Rc[X,])2] = E[(fm[A,])2] , ERc[X,]Im[X,] = 0. 

A number of important properties hold for these models, which we describe now. 
An excellent reference, with complete proofs is ,3J. 

We write U(z, r) for the open disk {w e C : \w — z\ < r}. The natura l do - 
main of convergence for f^{X,z) is U(0,/Ok), a.s., where p^ = |k|^^/^. By (2.f ), 
^[f^{X,z)f^{X,w)] = 0, and 



Uix,z) -- 


= ^an^^XnZ 


n 




n=0 


n 


[1-(J-1)'«1 
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1/2 



(2.2) ^UX,z)UX,w)] = Q^{z,w) := 



(1 + Kzw)^/" for K 7^ 0, 
e^^ for K = 0. 



The function Q^ possesses important symmetries. For |u| < p^, consider the Mobius 
transformation 

*"^^^ " 1 + KUz' 

which is a univalent mapping of U(0, p^.) to itself. This is an isometry relative to a 
metric with Gauss curvature 4k. Moreover, 



QAKi^), $"H) = A:{z)Aliw)Q^iz,w), 

where 

_ j {I + K\u\^y/'^^'^^1 + Kuz)-^/" for K 7^0, 
I exp(i|Mp — uz) for K = 0. 
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Figure 1. To study the zeroes in a neighborhood of a point u 
near the boundary of U(0,/9k), take the image under the map $", 
which maps u to 0, and consider the positions of the zeroes under 
this mapping. This is described in Remark |2.2[ 



While Qk is not invariant with respect to the isometrics $", one says it is covariant 
because of this property. Also note that for any z E 1LJ(0,Pk), 

\K{z)\^Pk as \u\^p^. 

Taking u to the boundary of the domain U(0, p^), conformally maps neighborhoods 
of to domains approaching the boundary. 

Gaussian analytic functions (GAF's) are important special cases of RAF's. Their 
zero sets have been studied in [7 , and many interesting questions about these zero 
sets continue to be studied. The reader may consult that reference and references 
therein. Our main result proves convergence in distribution of the zero sets of the 
RAF's, for a sequence of neighborhoods converging to the boundary. 



Theorem 2.1 (Main Result). Suppose that ATojATi,.. 
random variables with mean zero and satisfying (2.1) 



are i.i.d., complex-valued 
Let Z ~ (Zg, Zi, . . . ) he 
i.i.d., complex Gaussians with density TT~^e~'^' '^ on the complex plane. For each 
K < 0, and any continuous function Lp whose support is a compact subset o/U(0, p^), 
the random variables 

4:/,(X,C)=0 

converge in distribution, in the limit \u\ ^^ p,^, to the random variable 



E ^(0 ■ 



C:/„(Z,C)=0 

Since Qk is covariant with respect to the mappings <i>J^, and since the distribution 
of a Gaussian process is determined by its covariance, the zeros of the GAF, {^ : 
fK,{Z,(,) = 0} is a stationary point process with respect to these mappings. 

Remark 2.2. The mapping $^ was defined so that $"(u) — 0. Therefore, mapping 
the zeroes by $" maps the zeroes in a neighborhood of u to a neighborhood of 0. 
We have an illustration in Figure [l] to indicate this. The test function (/? is nonzero 
only in a window around 0. 

There is some interest in also considering random polynomials with real coeffi- 
cients. This motivates considering random real analytic functions. 
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Theorem 2.3. Suppose that Xq, Xi, . . . are i.i.d., real-valued random variables 
with mean zero and variance 1. Let Zq^ Zi, . . . he i.i.d., real-valued standard, normal 
random variables. For each n < Q, and any continuous function Lp whose support 
is a compact subset o/U(0,/9k), 



E 



v(*J're(0) converges in distribution to j, viO 



in the limit \r\ — > p,^, along any sequence satisfying the constraint r € M. 
3. Proof of the Main Result 

The proof is elementary. It uses several tools from probability theory and com- 
plex analysis such as the Lindeberg-Feller condition and Hurwitz's theorem. Since 
these are well-known to probabilists, these results were merely referred to in an 
implicit way in the version of our paper 5] in order to shorten the presentation. 
But here, we will also briefly review those tools. (In an earlier version of our paper, 
which is available in preprint form, we did also include these tools.) 

We transcribe the Lindeberg-Feller version of the central limit theorem from 
Durrett's textbook [2]. (See page 110 in the latest, online version.) 

Theorem 3.1 (Lindeberg Feller Theorem). Suppose that for each n, there is a 
sequence of real random variables (X„ ,„) such that for a fixed n these random 
variables are independent (for different m indices) and EX„ .„ — for all m. 
Suppose 

(i) Em=l ^Xl,m ^ CT^ > as n ^ C50, 

(ii) For all e > 0, lim„^oo E(|X„^,„|2. |x„ „,| > e) == 0. 

Then Sn — X]m=i ^n,m Converges in distribution to ax, where x is a standard 
normal (real) random variable. 

Durrett also lists an exercise in his textbook to derive "Lyapunov's theorem" 
from this. (It is Exercise 3.4.12 in the most recent version.) Let us also state this 
corollary. 

Corollary 3.2. Suppose that Xq,Xi, . . . are i.i.d., complex valued random vari- 



ables with mean zero, satisfying condition {2.1). Suppose that (a„.fc) is a complex 
valued set of numbers, for n G N and k G {0, 1, . . . }, satisfying 

i'^) \ SfeLo lo^n.fcP -> cr^ > as 71 —> oo, and 

(t') SfcLo |Q^n,fc|^ ^0 as 71 —> oo for some p>0. 
Then S'„ = X]fe=i RA'^n,kXk\ converges in distribution to ax- 

We will not prove these results, which are available in [5]. However, we will use 
Lyapunov's condition for the Lindeberg-Feller theorem to prove: 

Lemma 3.3. Suppose that Xq^Xi, . . . are i.i.d., complex-valued random variables 



with mean zero and satisfying (2.1). Let Z = (Zq, Zi, . . .) be i.i.d., complex Gaus- 
sians with density Tr~^e~'^' '^ on the complex plane. Then for any N G N, any 
zi, . . . , zjv G U(0, Pk.) and any Ai, . . . , Ajv G C, the random variables 

N 



Ea 

k=l 



f,iX,^:{z,)) 



converge in distribution, to the random variable J2k=i ^kfKiZ,Zk), as \u\ -> p^,, 
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Proof. We may write 

where we have left the dependence of a„.„(u) on Ai, . . . , A„ and zi, . . . , z„ hnplicit 
for the coefficients 

fc=i ^'^y^i') 
Since Q,^ is covariant with respect to the transformations $", this imphes that 

oo N 

n=0 i,k=l 

For the same reason the variance of the random variables in question is a constant 
function of u. 

To apply the Lindeberg-Feller conditions, wc need to check Lyapunov's condition 
for the sequence a„^fc(u). We note that for p = 4 

n=0 ■ n=0 ' "^^ '^'' 

Cauchy's integral formula implies 






for each fixed z e 1LJ(0, p^)- For k < the series sums to 



dC, 



1/k 



The second factor on the right hand side has norm bounded by 1 for all C, G C{0, 1). 
Moreover since |$^(2:)| converges to p^ — |k|~^/^ in the limit |m| — > p^, the second 
factor converges pointwise to in that limit, for every ( S C(0, 1) \ {1}. For k ~ 
the series sums to 

ei^i%xp((c-i)i$::(z)p). 

But for K = 0, we know po — oo and |$o(-2^)l — > oo in the limit |u| ^- oo. Since the 
real part of (C — 1) is non-positive, the same conclusion follows. In either case, the 
dominated convergence gives the desired result. □ 



Lemma 3.3 implies that the random analytic functions [A"(z)]^^/k(JC, $"(z)) 
converge in distribution to the random analytic function /„(.Z,z), in the limit 
1^1 — > Pk, in the sense that the finite dimensional marginals of the function values 
converge. This also implies convergence in distribution of the zero sets. We will 
write => for convergence in distribution. A clear and elegant proof of this fact 
has been provided by Valko and Virag in a recent paper they wrote on random 
Schrodinger operators [5]. 
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Lemma 3.4 (Valko and Virag, 2010). Let fn{^^ z) be a sequence of random analytic 
functions on a domain D (which is open, connected and simply connected) such 
that E/i(|/„(z)|) < g(z) for some increasing unbounded function h and a locally 
bounded function g. Assume that fn{z) => f{z) in the sense of finite dimensional 
distributions. Then f has a unique analytic version and fn ^ f in distribution 
with respect to local-uniform convergence. 

Because of this result we see that [A"(z)]~^/k(X,$^(z)) converges in distribu- 
tion to fK,{Z,z), with respect to the local uniform convergence. To converge in 
distribution with respect to local-uniform convergence means that for any function 
T oi f which is continuous with respect to the local-uniform topology, the random 
variables J^{fn) converge in distribution to J^{f). li Lp is & continuous function, 
compactly supported on U(0,Pk), then defining the function 

?:/(C)=0 

this is continuous with respect to the local-uniform topology. For the reader's 
convenience we paraphrase Hurwitz's theorem from Saks and Zygmund ^, page 
158. 

Theorem 3.5 (Hurwitz's theorem). If a sequence {hn{z)) of functions, continuous 
on a closed set K. and holomorphic in the interior of IC, is uniformly convergent on 
this set, and if the function h{z) = lim„_>.oo hn{z) vanishes nowhere on the boundary 
of the set IC, then, beginning from a certain value ofn, all the functions hn{z) have 
in the interior of IC the same number of roots as the function h{z) (counting every 
root as many times as its multiplicity indicates). 

From this theorem, it is easy to see that Af^p is a continuous function. This follows 
from the usual method of approximation by simple functions based on disks. 
Note that A"(z) is finite and non-vanishing for z e U(0,Pk)- Therefore 

Af^{[A-{z)]-^f^{X,<i>-{z))) = Af^{f^{X,^':{z))) . 

Moreover, /k(X, ^]t(z)) = means that z = $^"(i(;) for some w such that /k(X, w) = 
0. Therefore, we see that convergence of the distributions ol N^ {fi^{X, $"(z))) re- 
ally means that 

y^ '^'(^"(O) converges in distribution to \J tp{£^) , 

in the limit lul — > p^, as claimed. 



4. An Application 

Recently there has been some interest in looking at all random polynomials of a 
given degree with coefficients in {-hi, — !}• John Baez reported on this on his blog, 
and there was a popular article on the topic with numerical results obtained by 
Sam Derbyshire [T]. We can write such polynomials as 

n 

Pn{Xo,...,X„;z) =^Xkz'' , where Xq, . . . ,X„ e {+1, -1} . 

fc=0 
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13 



■'-2 -1.5 -1 -0.5 0.5 1 1.5 2 

(b)W8 



Figure 2. (a) This is the set of ah z e C, such that 
Pn{Xo, . . . , Xn', z) — for some choice of Xq, . . . , X„ e {+1, —1}, 
for n — 13. (b) The set of ah roots for polynomials with coefhcient 
Xq, . . . ,Xn E {l+i,l — i,—l + i,—l — i} ioT n = 8. (The roots are 
computed numerically using Octave's root program. Better figures 
are available on the arXiv version of this paper.) 



In Figure [2J we have plotted a representation of the zeros for 
sequence Xq, Xi,- ■ ■ e {+1, —1}, it is easy to see that 



Pn{Xo,...,Xn;z) -)■ /_i(X,z), 



13. For any fixed 
locally, uniformly 



as n — )■ oo. In particular the zero sets converge. 

Note that since the coefficients are strictly real, our main theorem does not apply. 
But this is why we stated Theorem |2.3| The proof of Theorem |2.3| is identical to 
that of Theorem |2.1[ The only difference is that the covariance changes. Instead 
of having E[f{X, z)f{X, w)] — 0, we have 



/(X,z) = /(X,z), 

because all coefficients are real. The Mobius transformations <I>J1 for r e M also pre- 
serve this property. That is why we restricted to those isometrics in the statement 
of Theorem [231 

One can also consider the polynomials p„(Xo, . . . , X„; z) where Xq, . . . , X„ are 
i.i.d., uniform in the set {1 + i,l — i, —1 + i, —1 — i}. In this case. Theorem 2.1 
applies. The simplest non-Gaussian RAF is given by 

oo 

f-i{X,z) = Y.X,,z\ 

n=0 

where Xq, Xi, . . . are i.i.d., random variable chosen from the set {1, —1, i, —i} with 
equal probabilities. Our theorem implies that the zero set of such a RAF near the 
unit circle is has a distribution which is close to that of the corresponding GAF. 
Let us define 



Zn = [j {ze 

,Xr,e{+i,-i} 



Pn{XQ,. ..,Xn;z) = 0}. 



A'n 



8 SHANNON STARR 

and 

Wn = U {zee : p„(Xo,...,X„;z) = 0}. 

Xo,...,x„e{i+'ia-'i,-i+i,-i-i} 

Baez notes several "holes" in the sets Z„ centered at points on the unit circle, such 
as 1 and — 1 . From Figure [2| one also sees several holes in the set W„ for n ~ 8 
along the unit circle. Due to our theorem, we can deduce that the holes in Z„ at 
±1 must close up. Similarly, all the holes in the set Wn along the unit circle must 
close in the limit n — > oo. 

For Zq,Zi,... i.i.d., real-valued standard, normal random variables, the zeros 
of /_i(.Z, z) are asymptotically uniformly distributed along the unit circle due to 
the result of Ibragimov and Zaporozhets, for example. The intensity measure of 
f-i{X,z) is asymptotically close to that of /_i(.Z, z), in the weak topology, near 
±1. Therefore, there cannot be a hole at 1 or —1. This is the simplest consequence 



of Theorem 2.3 Similarly, in the complex case. 

One can also consider this from another perspective. For a fixed value of z, we 
may define the set of function values for all possible coefficients 

Cn{z) = {pn{{xo,...,Xn),z) : a;o,...,a;„ G {1,-1}}. 

This satisfies the recurrence relation 

C„(z) = {1 + zw : w e C„_i(z)} U {-1 + zw : w e C„_i(z)} . 

Defining the set 

Ciz) - {/-i(y,z) : 2/o,yi,---e{i,-i}}, 

one can see that C{z) = {1 + zw : w E C{z)} U { — 1 + zw : w e C{z)}. This 
implies the Hausdorff dimension satisfies the bound 

log(l/|z|) 

For r e M, we know that C{r) C M. Therefore, it seems reasonable to conjecture 
that 

dimCfz) = fmin{2,log(2)/log(l/|z|)} for z € C \ K, 
^""^ [min{l,log(2)/log(l/|z|)} for z G M. 

It is easy to see that C(l/3) is the middle-thirds Cantor set; whereas, C(z/-\/2) 
is the rectangle {x + iy : a; G [—2,2], y G [— \/2, \/2]}. In Figure [s] (a) we have 
displayed Ci5(z) for z = e"^/V^/2. 

The fractal dimension of C{z) may pertain to Baez's conjectures that the zero 
sets {z : G Cn{z)} have a multi-fractal structure associated to z. In analogy to 
the dimension of varieties of smooth curves, one might guess that "typically" if z 
satisfies C{z) 3 then the set of w's in a neighborhood of z with C{w) B has 
dimension "approximately" equal to dim C(z). 

Similarly, defining 

B„iz) = {p„((xi,. ..,x„),z) : xo,Xi,- •• G {1,-1, i,-i}} , 
we have Bn{z) = Uue{i,-i.i,-i}{^ + zw : w e i?„_i(z)}. Defining 
B{z) = {f-i{x,z) : Xo,Xi,--- e {1,-1, i,-i}}. 
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Figure 3. (a) The set 6*15(2) for z = e"l^l^/2. (b) The set 
-88(2;) for z = e™l^\/2. (Better figures are available on the arXiv 
version of this paper.) 



this implies that <XvcaB{z) < 21og(2)/log(l/|z|). It seems reasonable to guess 
that dim. B{z) = min{2, 21og(2)/log(l/|2;|)}. An easy calculation shows -6(1/2) = 
{x + iy : x,y e [-2, 2]}. In Figure|] (b) we have plotted Bs{z) for z = e'''/^/2. 

One can also ask, for z such that B{z) 9 0, about the structure of the coefficients 
{{Xo,Xi,...) : Xo,Xi,-- • e {1,-1, 1, -i}} satisfying /_i(X,2;) = 0. One guess is 
that as \z\ -^ 1, the distribution converges in some sense to "uniform" with the the 
density appropriate for the intensity measure of the zeroes at that point. This guess 
is affirmed for GAF's. It is reasonable to conjecture that this is also true for RAF's 
with discrete coefficients. An interesting possibility is to relate the correlations of 
zeros for /_i(X, z) for a typical value of X, to the structure of the coefficients X 
such that f-i{X,z) = for a typical value of z. This may be a topic for further 
study. 
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